Surface magnetic canting in a ferromagnet 
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The surface magnetic canting (SMC) of a semi-infinite film with ferromagnetic exchange interac- 
tion and competing bulk and surface anisotropies is investigated via a nonlinear mapping formulation 
of mean-field theory previously developed by our group [L. Trallori et al, Int. J. Mod. Phys. B 
10, 1935-1988 (1996)], and extended to the case where an external magnetic field is applied to the 
system. When the field H\\ is parallel to the film plane, the condition for SMC is found to be the 
same as that recently reported by Popov and Pappas [Phys. Rev. B 64, 184401 (2001)]. The 
case of a field H± applied perpendicularly to the film plane is also investigated. In both cases, the 
zero-temperature equilibrium configuration is easily determined by our theoretical approach. 

PACS numbers: 75.70.Rf, 75.70.-i, 75.30.Kz 



Recently Popov and PappasS - motivated by their 
previous valuable experimental workn on the magnetic 
properties of an ultrathin (1.5 monolayers-thick) Fe film 
grown on the surface of a 15nm-thick Gd(0001) film - de- 
termined the zero-temperature phase diagram of a semi- 
infinite Heisenberg ferromagnet, with exchange constant 
J > 0, subject to a surface anisotropy, Ks, competitive 
with the bulk one, Kb- In the mean-field approxima- 
tion, since the system is inhomogeneous only along one 
direction (i.e., the normal to the surface), the problem 
is reduced to consider the energy of the one-dimensional 
modeltH 



E = -JJ2 cos (0« " °n+i) + K s sin 2 9 1 

71=1 

oo oo 

+ K B ^2 sin 2 0„ - 9HbH\\ ^ cos0„ (1) 



where 6 n is the angle formed by the classical vector mo- 
ment of the n-th layer with the film plane. To fix ideas, 
in the following we will assume Kb > (i.e., the bulk 
anisotropy favours the alignment of the vector moments 
along the film plane, 6 — 0) and Ks < (i.e., the surface 
anisotropy favours the alignment along the film normal, 
9 = tt/2). Hn denotes an external magnetic field applied 
parallel to the film plane. 

Expanding the energy in Eq. (|l|) to second order for 
each n , E Eq + 9 T A 0,_and performing a stability 
analysis, Popov and Pappaso showed that it is possible 
to have surface magnetic canting (SMC) provided that 



1 + h\\ < 



(ks ~ 1) 



(2) 



where ks — 2Ks/J, ks = 2Kb /J and h\\ = gfiBH^/ J. 
Hence it follows that for ks < — 1 — hn, the surface is 
always canted whatever the value of the in-plane bulk 
anisotropy; otherwise there is SMC for \Ks\ exceeding a 
threshold value which depends on Kb, J and Hn. 



Surely, the theoretical analysis of Popov and PappasS 
is quite original and valuable. However, it is worth not- 
ing that, in the special case H\ = 0, a condition for sur- 
face magnetic canting analogous to Eq. was obtained 
by our group some years agoO using a rather different 
method: i. e n( -a-|H|©iiiinear mapping formulation of mean- 
field theory.QBBfl&El Within this framework, the proper- 
ties of a magnetic film are formulated as a problem in 
nonlinear dynamics, in terms of an area-preserving map, 
where the surfaces are introduced as appropriate bound- 
ary conditions. In this Brief Report, using the method 
developed in Refs. |,|§,@,1,|, not only we obtain the 
condition for surface magnetic canting, but also we cal- 
culate the zero-temperature magnetization profile of the 
film, both in the caser-pf the model, Eq. (|l|), considered 
by Popov and PappasEJ, and in the case that a magnetic 
field is applied perpendicularly to the film plane, see Eq. 
@ belowB 

Let us summarize our method. The equilibrium con- 
figurations of the semi- infinite [n > 1) film are obtained 
from Eq. (1) by ^-derivation 

0E 

J sin(6»„ - n+ i) - J sin(0 n _i - 0„)(1 - S n ,i) 



d9 n 

+ Kssm(29 n )S nil +K B sm(29 n )(l- 
+ g(J.BH\\ sm9 n = 

Introducing the variables s„ = sin(0 n — 6 n —i) 
(n > 2) we obtain the nonlinear mapping 

Js n+ i = Js n + K B sin(29 n ) + gn B H\\ 

n +i = On + sin _1 (.s„ + i) 

while on the surface plane n = 1 we have 



(3) 

in the bulk 



sin 9 r , 



(4) 



Js 2 = 
9 2 = 



The map, Eq. (g) 
determinant is 1 

detJ = dct 



-Kssin(20i) + g^BH\\ sin0i 
1 +sin- 1 (s 2 ) (5) 

is area-preserving because its Jacobian 



,9s „ 



»0„ 
ftW, H . 



2 



det 



S(6 n ) 
1 



(6) 



where <S(0 n ) = ^ cos(26»„ 



^^cos(9„. The trajec- 



tories in (0, s) space are associated with equilibrium con- 
figurations, while the fixed points of the map correspond 
to uniform ground states of the bulk system. For K$ < 
and Kb > 0, the hyperbolic fixed point: (0oo, Soo)=(0, 0) 
represents the energetically stable bulk configuration, 
while the elliptic fixed point: (9oo, s <x> )=(n/2, 0) repre- 
sents an unstable bulk configuration. 

Now we observe that, by introducing a fictitious plane 
n = 0, characterized by the angle 9q and the variable 
S\ = sin(#i — 9q), the nonlinear mapping Eq. (Q), valid in 
the bulk (n > 2), can be assumed to hold even for n = 1, 

Js2 = Js\ + Kg sin(2#i) + g[igH\\ sin#i 
9 2 = 6 1 +sm- 1 (s 2 ) (7) 

provided that the following boundary condition is satis- 
fied 



Ks - Kb . , na , 
Si = sm(26q). 



(8) 



Thus, among all trajectories in (6, s) space obtained from 
Eqs. (H)(0), only those which satisfy Eq. (||) represent 
equilibrium configurations for the semi-infinite ferromag- 
net. 

By linearizing the map in the neighborhood of a fixed 
point JlJ one is led to solve the eigenvalue equation 



A 2 - Tr J A + det J = 



(9) 



Near the hyperbolic fixed point (9oo, Soc)=(0, 0), one ob- 
tains two real eigenvalues (Ai < 1, Aa > 1) 



Al.2 =1 



S(8 C 



(10) 



and the slopes, in the (0, s) phase space, of the orbit in- 
flowing to (" 1" ) and outflowing from (" 2" ) the hyperbolic 
fixed point are, respectively 



mi.2 



A, 



(5s „ 



(11) 



(0,0) 



Within this theoretical framework, the condition for 
surface magnetic canting is that m s , the (negative) slope 
of the boundary condition curve Eq. @, calculated in 
the hyperbolic fixed point 



(0,0) 



2(K S - K B ) 
J 



(12) 



should be smaller than m 1; the (negative) slope of the 
trajectory inflowing to the hyperbolic fixed point, see 



FIG. 1: Map phase portraits calculated from Eqs. (|4|)(FJ us- 
ing Hamiltonian parameters J — 1, Kb = 0.1, Ks = —0.3 
for different values of a magnetic field H« applied parallel to 
the film plane: (a) g^sH^ — 0; (b) 0.05; (c) 0.155 (threshold 
value, see Eq. ([lq); (d) 0.3. The dashed curve represents the 
boundary condition at the surface, Eq. (^). Arrows denote in- 
flowing and outflowing trajectories associated with hyperbolic 
fixed points. 



Eq. (|l l|) . Thus, the condition for surface magnetic cant- 
ing turns out to be 



2(K S - K B ) 
J 



< 



S(#oo) 



<S(0 a 



(13) 



Taking into account that Ks < and Kb > 0, this 
equation can be rewritten as 



\K S \ > 2K B +gn B H\\ 



A.J 



4J 

2Kb - 9HbH\\ 



2Kb + 9^bH« 



4J 



(14) 



which is readily seen to be completely equivalent to 
Eq. @. Thus, in the case of external magnetic field 
parallel to the film plane, for sufficiently high surface 
anisotropy (i.e., 2|ifg| > J + gfigHu), the surface is al- 
ways canted whatever the value of Kb ■ Otherwise, there 
is SMC for \Ks\ exceeding a threshold value which de- 
pends on Kb, J and Hn: see Eq. (|l4|). 

In Fig. 1 we report the different map phase portraits 
obtained from Eqs. (§)(0) for different, increasing val- 
ues of H\\ . The boundary condition curve - which does 
not depend on the value of Hu, see Eq. (g) - is also re- 
ported (dashed line). In zero field, H\\ — (see Fig. la), 
for the chosen values of the parameters, the condition, 
Eq. (|l3|), for the slopes m s and mi calculated in the hy- 
perbolic fixed point, is satisfied: thus, the configuration 
with surface magnetic canting is the ground state. As Hn 
is increased, the slope mi of the trajectory inflowing to 
(0oo, Soo) =(0,0) becomes more negative and 9\ decreases 
(see Fig. lb). Finally, for H\\ > Hc\\> where Hq\\ is a 
threshold value given by 



9^bH c \\ = 



[2(K B - K s )} 2 
J + 2{K B - K s ) 



2K f 



(15) 



the boundary condition curve is no more intersected: 
thus, the uniform configuration with all spins parallel 
to the film plane becomes energetically favoured (see 
Figs. lc,d). 

Let us now consider a semi-infinite ferromagnetic film 
(J > 0) with competing surface and bulk anisotropics 
(Kg < 0, Kb > 0) and with an external magnetic field 



applied perpendicularly to the film plane.ti] The energy 



0" 



3 



FIG. 2: Map phase portraits calculated from Eq. ([17]) using 
Hamiltonian parameters J = 1, Kb = 0.1, ifs = —0.3 for dif- 
ferent values of a magnetic field H± applied perpendicularly 
to the film plane: (a) gfi B H± = 0; (b) 0.10; (c) 0.20 (thresh- 
old value, g/J-BHci- = 2Kb)', (d) 0.3. The dashed curve repre- 
sents the boundary condition at the surface, Eq. (|^). Arrows 
denote inflowing and outflowing trajectories associated with 
hyperbolic fixed points. 



is 

oo 

E = - jJ2 cos (0« - e n+i) + K s sin 2 9 1 

n=l 

oo oo 

+ Kb^sw 2 6 n - g/j, B H ± *£2 sin 6 n (16) 

n— 2 n— 1 

The map equations are now (n > 1) 

Js n+ i = Js n + K B sin(20„) - g/J, B H ± cos 6 n 
n +i = On + sin _1 (s n+ i) (17) 

while the boundary condition at the surface plane turns 
out to be the same as in the case of in-plane magnetic 
field, Eq. (||). The Jacobian determinant takes the same 
form as in Eq. (^): the only difference is that now S(6 n ) — 
^fi-cos(20„) + smH± sm g n , j n contrast with the case 
of in-plane field, now the hyperbolic fixed point (9^, s^) 
is characterized by an angle, #oo, which depends on the 
field intensity 



gp, B H± 



2K E 



, for gfi B H± < 2K B , 
, for g/J, B H± > 2K B 



(18) 



and Soo — 0. In Fig. 2 we report different map phase 
portraits, obtained from Eq. (|L7]) calculated for different, 
increasing values of H_\_. The boundary condition curve 
Eq. (||), independent of the value of H±, is also reported 
(dashed line). It is apparent from the map topology (see 
Figs. 2a,b) that for < H± < Hc±, where g/i B Hc± = 
2K B is a threshold field value, one does have surface 
magnetic canting even in the case that SMC is absent for 
zero field. In contrast, for H± > Hq± (see Figs. 2c, d), no 
intersection is possible between the boundary condition 
curve and the trajectory inflowing in the hyperbolic fixed 
point (6oo, Soo)=(f , 0): thus, the uniform configuration 
with all spins perpendicular to the film plane becomes 
the ground state. 

In Fig. 3a,b we report the zero-temperature equilib- 
rium configurations of the semi-infinite film obtained via 
the map method for different values of the intensity of 
an external magnetic field, applied in-plane or perpen- 
dicularly to the plane, respectively. In both cases we 
have chosen Hamiltonian parameters which, in zero field, 
give rise to surface magnetic canting. We observe that 
as the in-plane field H\\ is increased (see Fig. 3a from 
top to bottom), the surface angle d\ decreases until, for 
H\\ > Hc\\ i a uniform configuration with all spin in-plane 



FIG. 3: (a) Equilibrium configurations of the semi-infinite 
film with Hamiltonian parameters J = 1, Kb ~ 0.1, Ks = 
—0.3, as calculated from the map equations Eqs. (^)(§) and 
boundary condition at the surface Eq. (g|). The different 
curves refer to selected values of g^BH^ — 0, 0.10, 0.15, and 
H\\ > Hc\\, see Eq. (|l5|), with the field increasing from top 
to bottom. The lines are guides to the eye. (b) The same 
as in (a), but for magnetic field applied perpendicularly to 
the film plane. The different curves refer to selected values 
of g/i B H ± = 0, 0.10, 0.15, 0.19, and H± > H c ± (where 
gfiBHc± = 2K B = 0.2), with the field increasing from bot- 
tom to top. 



is obtained. In contrast, in the case of increasing perpen- 
dicular field (see Fig. 3b from bottom to top) , the surface 
angle 6\ increases and so does the bulk angle 0^, defined 
As gfi B H ± 



in Eq. (18). As gfi B H± reaches the threshold value 2K B , 
a uniform configurations with all spins perpendicular to 
the film plane is obtained. 

It is worth noting that by the map method it is possible 
to calculate, with any desired accuracy, the equilibrium 
configuration even in the case of a film with a finite num- 
ber TV of spins, provided that the boundary condition on 
the second surfacejju sjv = — Ks ~ K j?_ sin(26*Af), is taken 



into account, in addition to Eq. (|j).El 

Finally, we observe that the approach described here 
can be used to investigate the zere-temperature mag- 
netic properties of the Fe/Gd film,H provided that the 
model (0) is modified in order to account for the anti- 
ferromagnetic couplingJaetween the Fe surface plane and 
the Gd underlying oneJ13 However, owing to the absence 
of frustration effects, only qualitative modifications to 
the results obtained here are expected. Moreover, it is 
worth noting that the non-uniform ground state is not 
affected by a surface enhancement of the exchange con- 
stant (which is present in the Fe/Gd film, since the 1.5 
monolayers of Fe order at a Curie temperature substan- 
tially higher than Gd). In fact, since at T — the mag- 
netization of each plane takes the saturation value, the 
surface exchange only keeps the spins on the surface par- 
allel to each other, without modifying their orientation 
with respect to the spins on the underlying plane: i.e., at 
T = the possibility of a non-uniform ground state only 
arises from the competition between surface and bulk 
anisotropics In contrast, the interpretation of the quite 
interesting finite temperature properties of Fe/Gd films, 
like the two-step reorientatipn transition experimentally 
observed by Arnold et al.n would require not only an 
improvement of model (|l|) but also a much more refined 
analysis. In fact, at T ^ both the modulus and the ori- 
entation of the magnetization vary with the plane index, 
so that the dimensionality of the map increases and the 
equilibrium configuration of the film cannot generally, be 
calculated without resorting to some approximation.Q 

In conclusion, a mean-field theoretical method, where 
the equilibrium properties of a magnetic film are formu- 
lated in terms of an area-preserving map and the sur- 



4 



faces are introduced as appropriate boundary conditions, 
was exploited in this paper. In particular, the model of 
a semi-infinite Heisenberg ferromagnet with competing 
bulk and surface anisotropics was investigated when an 
external magnetic field is applied parallel or perpendic- 
ular to the film surface. In both cases, the condition for 
surface magnetic canting as well as the zero-temperature 
equilibrium configuration were easily determined by such 



a theoretical approach. 
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